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ABSTRACT 

During the process of collapse of a massive star, a cavity is generated between 
the central iron core and an outer stellar envelope. The dynamics of this cavity, 
filled with plasma and magnetic field of the rapidly rotating proto-magnetar's 
magnetosphere, is believed to be very relevant in understanding supernovae and 
gamma-ray bursts. The interactions of the pressurized conducting plasma and the 
magnetic fields are described by a set of magnetohydrodynamic (MHD) equations 
with poloidal and toroidal plasma flows not aligned with magnetic fields. A 
sequence of MHD equilibria in response to the increasing plasma pressure in 
the cavity, by continuous filling from the rotating magnetosphere, is solved to 
account for asymmetric supernovae, highly collimated gamma-ray burst jets, and 
also active galactic nucleus plasma torus. It is shown that the magnetosphere of 
the central compact star is likely the central engine of supernova and gamma-ray 
burst by feeding them plasma, magnetic energy, and rotational energy. 

Subject headings: supernova, gamma-ray burst, magnetic towering 
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Introduction 



Ever since Baade and Zwicky had coined the term super-novae and had suggested super 
nova process as the end of an ordinary star to become a neutron star (IBaade fc Zwicky 
19341 ). supernovae have been puzzling scientific minds for decades. They are believed to be 
the catastrophic end of massive stars' life cycle. Supernovae other than the type la standard 
candle, where the entire st ar is incinerated, are belieyed to have gone through a gravitational 



1990; iKotake et al. 



20061 ). The progenitor iron core with 



collapse of the iron core (jBethd 
about 10 4 Km radius is enclosed by the stellar hydrogen envelope extending out to some 
10 8 Km radius. The collapse of this progenitor core is governed by the sound speed profile, 
which decreases as the radial position increases. The infalling velocity in the outer part of 
the core is supersonic with respect to the local sound speed, while the inner part is subsonic 
with the interface at about 300 Km radius. Because of this condition, the outer part tends 
to pile up while the inner part continues to free-fall. This generates a proto-neutron star 
plus the overlying stellar material. A cavity is formed between the proto-neutron star and 



jet formations was investigated i 
associated with nuclear detonati 


LeBlanc & Wilsor 


1970 


), and an e 


on was proposed ( 


Bisnovatvi 


1971) 



Some of the core-collapse supernovae appear to be associated to gamma-ray bursts 



( Galama et al. 



1999 



Campana et al.ll2006 



Pian et al.ll2006 



Soderberg et al. 



2006 



Mazzali et al 



20061 ). As a result, there could be a very close relationship between supernovae and gamma- 
ray bursts. By the extraordinary energy output, servera l scenarios have been contemplated 
as the progenitors of gamma-ray bursts ( iMeszarosI 120021 ) . For massive progenitor stars with 
more than 14 solar masses, like rapidly spinning Wolf-Rayet stars, the inner core could 
promptly collapse to a black hole circumscribed by a massive torus, failing to generate a 
core rebound in the cavity. Accretion of the surrounding massive torus of nuclear density 
material at a later time drive an outburst along the rota tional axis breaking out the s tellar 
envelope. This is the Collapsar (or Hypernova) scheme ( iMacFadyen fc Woosleylll999l ). In- 
stead of promptly collapsing to a black hole, there also could be a two-stage collapse by first 
forming a rapidly rotating neutron star temporarily stabilized by rotation, an d later collapses 
t o a b lack hole after losing some angular momentum, which is the Supranova ( jVietri fc Stella 



19991 ) scheme for long gamma-ray bursts. When the inner core does not collapse to a black 
hole, it rebounces as the inner core reaches nuclear densities, generating a rebounding shock 
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in the cavity. The release of the gravitational binding energy through neutrin o bursts has 



been considered as the primary candidate in fueling the supernova explosion ( iBethd 11990 



Kotake et al. 



2006 



Bethe fc Wilson 



19851 ) . Nevertheless, as this outgoing shock meets the 



infalling outer very thick envelope, energy of the shock wave gets dissipated, and the shock 
is stalled. By which energy source and mechanism that the shock could be reignited is still 
an unsettled issue. 

Recently, polarimetry observations of supernova optical emissions have revealled dif- 
ferent degrees of polarization along a fixed axis of the supernova. In a very collisional 
environment of a supernova, this polarization result implies a nonzero volume average of 
the microscopically random electric field vecto r of each emission, indicating supernovae, or 



some them, are aspherical ( jHowell et al. 



2001 



Wang et al.l 120031 ) . In view of other astro- 



physical ejection events known to be powered by magnetic fields, these polarization ob- 



servations have given grounds t o reexamine the magnetic 



supernova (iWheeler et al. 



Burrows et al. 



2007 



200 



0. 



2002 



Arde 



Komissarov fc Barkov 



ian et al. 



2005 



iield as the central engine for 



Uzdensky &: MacFadyen 



2007 



20071 ). The primary concern in this renewed 



magnetic approach is the magnetic collimation mechanism. It has been proposed the pres- 
ence of a n accretion d i sk wi thin the c avity of the collapsing iro n core, by which magnetic 
towering ( Lynden-Belll 2003 ) and jets ( Blandford fc Payne 1982 ) could be generated along 
the rotational axis. 

Should magnetic field be the energy source of supernova, a simple estimation indi- 
cates that the surface fields of the rapidly rotating proto-neutron star would be around 
10 15 Gs, which qualifies it as a magnetar. Such magnetar scenario could be acomplished if 
the rotation of the proto-neutron star is fast enough. With millisecond periods, dynamo 



effec ts inside the pulsar could 



10 3 (iDuncan &: Thompson 



3 e triggered, enhancing the magnetic fields by a factor of 



1992 



Thompson fc Duncan! Il993l ). These magnetic fields would 



be launched from the magnetosphere together with the plasma to fill the cavity, leaving a 
normal spun-down neutron star at the center after supernova explosion. In general, the dy- 
namics of this ca yity driven by the rapidly rotat ing inner core is referred to as the 'Pulsar in 
a Cavity' model (jUzdensky fc MacFadvenl 120071). Simulation has been the principle tool i n 



investigating the dynamics in this cavity (iBurrows et al 



2007 



Komissarov fc Barkov 



2003). 



Most frequently, an accretion disk scenario with mag netic fields anchore d on it is considered. 
Under given initial conditions, magnetic towering (iLynden-Belll 120031 ) due to the angular 
momentum of the disk generates a collimating magnetic column along the polar axis. 
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Here, we seek to describe the cavity structures, filled with plasma and magnetic fields, 
through a sequence of axisymmetric magnetohydrodynamic (MHD) steady states (equilib- 
ria) with poloidal and toroidal plasma rotations as a response to the increasing plasma 
pressure. This type of steady state analysis has been used to gain important insights of dy- 
namic processes such as a strophysical jets in terms of spatially self-similar MHD equilibria 
( Blandford fc Payne 19821) and mag netic towering in terms of magnetostatic analysis of disk 



driven equilibria ( iLynden-Bell 



200 3J). The present problem differs from the relativistic pulsar 



wind problem wi th hot plasma (ILoyelace et al. 
similar analysis 



(jOkamoto 



1978 



Jprenderg 



ast 



Begelman 



2005 



1994 



1986 



Camenzind 



19861) especially under self- 



Gourgouliatos &: Vlahakisl |2010| ) and with cold plasma 



Okamoto 



20021 ) in that the spatial domain is bounded and 



very finite. 

In Sec. 2, the axisymmetric divergence-free rotational MHD formulation is briefly pre- 
sented, and in Sec. 3, by assigning two source functions, the Grad-Shafranov equation for 
rotational equilibrium is obtained. The nonlinear poloidal flux function is solved in Sec. 4 for 
asymmetric supernova configuration in the cavity, whereas cusp-like funnel polar collimated 
gamma-ray burst jet configuration is solved in Sec. 5. Active galactic nucleus plasma torus 
and some conclusions are finalized in Sec. 6. 
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2. Divergence-Free Rotational MHD 



The standard steady state MHD equations are 



V • (pv) 
p(v ■ V)v 

v x B 
Vxl 
V -B 



, 

t ft ^ GM - 

J x B — yp — p — — r 
-E , 
pf , 
, 
. 



Here, p is the mass density, v is the bulk velocity, J is the current density, B is the magnetic 
field, p is the plasma pressure, M is the mass of the central body, p is the free space 
permeability, and 7 is the polytropic index. To describe plasma equilibrium with poloidal 
and toroidal flows, this set of equations is quite inconvenient. We will use the MHD equations 
for divergence-free plasma flows, where the plasma velocity is density weighed so that the 
plasma density a ppears only through the density weig hed velocity w* = (ppY^v, except in 



the gravity term (ITsui et al. 



2011 



Tsui fc Navia 



20121 ). 



{pp) 1/2 V ■w, + tv,-V{pp) 1 ' 2 
x V x «J„ - B x V x B 

w*x B 
VxB 
V-B 
V • v 



GM. 



, 

V/ip* - p : 

T" 

{pp) 1/2 V§ = v$ 
pj , 
, 
, 



(2) 

(3) 
(4) 
(5) 
(6) 



where pp* = pp + w 2 /2 is the total plasma pressure. With divergence- free flows of Eq. El 
density flux conservation gives v- Vp = 0. By using the (w*, pp*) representation, this becomes 
w* ■ V(zip) 1 / 2 = 0, and Eq. [T]then becomes V • = 0. Furthermore, making use of Eq. HJ 
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the w* term and the B term are in symmetric form in Eq. [2j Because of Eq. the magnetic 
field can be represented through a vector potential. Under axisymmetry, this vector potential 
allows the magnetic field be represented by two scalar functions, which reads 



An ( 19* \ 

rsmt' \ r 06 or J 

j A ( + ^dF _8F _d 2 ^ _ 1 d 1 d$ \ 

rsin^ \ r 86 ' dr ' dr 2 r 2 86 sin 6 06 J 

Here, A carries the physical dimension of poloidal magnetic flux such that * is a dimen- 
sionless poloidal flux function, and A F is a measure of the axial plasma current. Likewise, 
because V • u>* = by Eq. [TJ we also have 



= A(,(W x V0 + F'V0) . (8) 

We note that An and A' are not the maximum values of poloidal magnetic and density 
weighed velocity fluxes. They are reference vlues only as such that the dimensionless poloidal 
flux functions are not normalized to unity. The condition ■ Vp = gives p = p(^f'). The 
scalar product of on Eq. [3] results in $ = $(*'). Consequently, with p = p(\& / ) , the right 
side of Eq. [3J can be written as the gradient of $ = $(*') only giving the second equality. 
Taking the scalar product again of B and to* on Eq. [3J give $ = $(*) and $ = $(*') 
respectively thus *' = *'(*). Substituting Eq. [7] and Eq. [8] to Eq. [3j gives 

(V*' x V$) x (V* x V$) - FV$ x (V*' x V$) + F'V$ x (V* x V$) 

= Iiv# . 

4,4, 

Taking the scalar product of V* on this equation gives 



-^--^W-V$ = (V0) 2 W • (F'W - FW') , (9a) 



7 



3. Rotational Grad-Shafranov Equation 



As for Eq. [5J with m' = m'(m), the <fi component gives 



l2 dm' ( dmdF' dmdF'\ 2 / dm dF dm dF\ 

~dm\~d7~dT + ~dl~dVJ^ \~d^~dJ + ~bJ~dV)~ ' 

We use Eq. [9b] to eliminate F' . It is noted that F and F' are not functions of m only for 
rotational equilibrium due to the (V0) 2 factor. With (V</>) 2 = (l/rsin^) 2 , we then get 



A' 2 (^—\ 2 - A 2 
<9\I> 



5$ dm <9F 
<9r <9# <9# dr 



A' dm' d$ 



d d 2 . 2 



By inspection, we have 



dm' 
K dm } 



A 



0n2 



A> 



F(r,6) = k{arf sin 2 # 
1 1 1 dm' 3$ 

+ ————— = 



A;a 2 A dm dm 



(10a) 
(10b) 



We note that a is a normalizing factor of r, and k is an inverse scale length such that the scalar 
function F(r,6) has the correct dimension of 1/r in Eq. [71 An arbitrary constant could be 
added to Eq. UOal However, this constant would correspond to an externally applied current 
along the axial direction, which would be appropriate for laboratory tokamak plasmas but 
not for astrophysical plasmas. Substituting d§/dm of Eq. UObl to Eq. l9bl gives 



F'M) = §F M ) + ^^M) = (J) . (ID 

Choosing the following linear source function leads to 



m'(m) = bm , 



(12a) 
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F'(F) 



b U'n 



aF 



(12b) 



We note that b is an independent model parameter, but a, denned in Eq. Il2bl is determined 
by b together with the parameters (Aq, A' q ). 

In order to analyse the 9 and r components, we choose to write the generalized total 
pressure in separable form 



HP*(r,9) = /ip*(r,*) = HPoP*i(r)p*2(^) 



(13) 



where p has the physical dimension of pressure p*, and p*i and p^ are dimensionless func- 
tions of r and \I> not normalized to unity, and p* is defined by 



HP* 
ka 2 



HP* 
.ka 2 



+ (At - A' W)F 



The 9 and r components then become 



1 



1 



(14) 



dr 2 



+ sin# 



5 



1 8^ 



89 V sin 89 



F 8 



((bA' ) 2 - A 2 ) ka 2 8$> 



(hp*) 



(15) 



F 8 



{{bA' ) 2 - A 2 ) ka 2 dr 



8^ 

dr 

(hp*) + 



dr 2 



8 



+ sin#— 



1 8^ 



89 \sin^ 89 
F 8^ 8 



(hp*) 



((bA' ) 2 - A 2 ) ka 2 dr 8V 

r sm 9jj,p — — 



((bA' ) 2 - A 2 



(16) 



The right side of the r component correspond to the explicit derivative 8/8r, the implicit 
derivative [8^ / 8r){8 / 8^) of generalized total plasma pressure, and the gravitational term. 
Making use of the 9 component to eliminate the implicit derivative term, the r component 
becomes 
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/W^P* 2 (*) +M*)^- = , (17) 



where we have defined a normalized radial coordinate z = ar. To satisfy this equation with 
p = p(\& ) , we take p*i(z) as below to get 



P*i(z) = \ , (18) 
P = ^M*) = PoP* 2 (*) = PoP* 2 W • (19) 

This shows that p is a funtion of \I/ only, consistent to the divergence-free flows, with po as 
the amplitude and p* 2 (\l/) as the dimensionless functional dependence. 

To solve the 6 component, let us further specify the second source function p* 2 as 



p„ 2 (tf) = (C± |^| 2m ) > . (20) 

Since plasma pressure is positive definite, \l/ 2m is tken in absolute value for any m. The 9 
component then reads 



2 + sin C -^75J im^sin 2 ^!^)!*! 2 ™- 1 = , (21) 
do \ smv ad J 



r 



dr 2 



n 2 PP0 

Pp 



{Al - (bA' )*)a* ■ 

This is the rotational Grad-Shafranov equation for divergence-free rotational flows with /3 P 
as the poloidal plasma (3. We remark that the static case can be recovered by taking A' Q = 
and with the corresponding limit of pp* and (3 P . This equation in spherical coordinates 
describes the structures in the cavity resulting from the interactions between the conducting 
plasma and magnetic field s. The correspondin g equation without rotational flows has been 



derived by many authors (ILovelace et al.lll98q ) in cylindrical coordinates 
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4. Asymmetric Supernovae 

Writing ty(r,9) = R(r)Q(9) in separable form, the rotational Grad-Shafranov equation 
becomes 



4^ + h^Te ) ^M^ier- = o . (22) 

We will arrange the above equation in the following form to solve the poloidal magnetic flux 
function ^ by separation of variables 



r>~ = -I S in S -K^f) Tm ^si„^p-, l( ,)^|e|— 

= n(n + l) . (23) 

The R equation gives R(z) — l/z n and R(z) = z n+1 as two independent solutions. As for 
the 9 part, in order to be separable, we choose R(z) — l/z n and p*i(z) — 1/z to get 



Sin ^f-^^V^ + 1 ) = Tmf3 p z" ( -) 2m 2 sin 2 9\Q\ 2m - 1 
d9 \sm9 d9 J K ' v \z n J 1 1 

= Tm{n) (3 P sin 2 9\Q\ 2m ^~ 1 , 
with 2n{m — 1) = 3 or m{n) = (3 + 2n)/2n. Defining x = cos 9, we have 



(1 _ X 2 )^M. + n (n + l)e(x) = T"i(n)/3 P (1 - x 2 )^! 2 "^- 1 . (24) 

This nonlinear equation with (3 P 7^ usually would have asymmetric solutions of 6(x). 
Nevertheless, we note that this equation is symmetric with ±x, hence symmetric nonlinear 
solutions of Q(x) are also allowed. As for the magnetic field components, they are now given 
by 



B r = -7T- , 25a 

z z ox 
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A a 2 



z (1 — x 2 ) 1 / 2 dz 



A a 2 



z (1 - x 2 ) 1 ' 2 a 



(25b) 
(25c) 



These magnetic fields are expressed in terms of the reference field a 2 A$. We remark that the 
Bq and B$ fields are scaled by the (1 — x 2 ) 1//2 factor in the denominator, which is singular 
at x = ±1 or 9 = 0, it. Since this singularity is quadratically not integrable, the magnetic 
energy would diverge at this location. To eliminate this divergence, we require the oscillating 
function Q(x) be null at the poles x = ±1. As a result, the oscillating Be and B^ fields grow 
to large amplitude near the poles, before they plunge to zero at the poles. The constraint 
that 0(±1) = makes /3 the eigenvalue of Eq. [24j The magnetic field lines and, in particular, 
the poloidal field lines are described by 



B r B e B^ 
dr rdO r sin i 
V(z,x) = R(z)e(x) = C 



(26a) 
(26b) 



where the poloidal field lines are given by the contours of ^(z, x) on the (r — 9) plane, which 
is also shared by the poloidal velocity stream lines. 

To solve Eq. [24j we note that, for (3 p = 0, G(x) is given in terms of Legendre polynomials 
by Eq. I29bl in Sec. 4 with n axisymmetric lobes, as shown in Fig.l with n — 1 and n = 3. 
However, (3 p = would imply = 0, or equivalently p = 0, w* = 0, and A = by Eq. [TU 
which amounts to no pressure, no flow, and no magnetic fields. For j3 p ^ 0, Eq. [241 defines the 
nonlinear Q(x) for a given n. To solve for equilibrium 0(x), we take the lower sign in Eq. [201 
and Eq. [24]with p^i^) — (C — \^/\ 2m ) to have a balance with the poloidal magnetic pressure. 
We start from x — — 1 with 0(— 1) = and a slope 0'(— 1) as such to reach x — +1 with 
G(+l) = 0. In the parameter space of (B'(— 1),/3 P ) for a given n, there will be (n-l)-lobe, 
(n-2)-lobe, 1-lobe structures in general. Rather than covering the (6'(— 1),/3 P ) parameter 
space, we take the derivative of the linear (3 P = solution of Eq. I29bl as the reference slopes 
at x — — 1, with 
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dP 

e'(x,n) = -2x^ . (27) 
dx 

For n = 2, a 1-lobe structure is shown in Fig. 2 with ©'(— 1, 1) and its corresponding (3 p = 9.57, 
where the lobe asymmetry is not evident. For n = 3, Fig. 3 shows a 2-lobe structure with 
0'(— 1, 3) and /3 p = 6.69, and for n = 3 in Fig. 4 is another 2-lobe structure with a different 
boundary derivative 0'(— 1, 2) and at a different pressure (3 P = 6.12. Naturally, there are also 
1-lobe solutions. For n = 4 with ©'(— 1,4), Fig. 5 shows a 3-lobe structure with j3 p = 6.05, 
which is rather symmetric as in Fig. 2. However, the symmetric lobes of Fig. 2 and Fig. 5 
are nonlinear lobes as they differ from the linear j3 p = corresponding lobes of Fig.l, 
where the (x,y) scales in the three figures have the same proportions. For n = 4 and 
under the same boundary derivative ©'(— 1,4), Fig. 6 shows a 2-lobe structure as pressure 
increases to (3 P = 6.41. In this figure, one lobe is overwhelmingly dominating over the 
other, as in Fig. 3 and Fig. 4, showing the nonlinear nature of the equation. For n = 4 
with 0'(— 1,3), Fig. 7 shows an asymmetric 3-lobe structure with (3 P = 12.3. For n = 4 
and changing the boundary derivative to ©'(— 1, 1), Fig. 8 shows another asymmetric 3-lobe 
structure with f3 p = 47.0. Comparing the 3-lobe structures of n = 4 in Fig. 5, Fig. 7, and Fig. 8, 
all with proportional (x,y) scales, the lobe amplitude decreases as the boundary derivatives 
decreases from ©'(— 1,4) to ©'(— 1,1). Multipole magnetic struct ures similar to ours have 



been explored in magneto- rotational supernovae ( jBisnovatyil 120081 ) 



We note that with R(z) = l/z n the poloidal magnetic flux function is ^(z, x) = 
R(z)Q(x). The poloidal magnetic field lines and the density weighed velocity stream lines 
with ^(z, x) = C are shown in Fig. 9 for n = 3, (3 P = 6.12, and in Fig.10 for n = 4, (3 P = 6.41, 
corresponding to Fig. 4 and Fig. 6 respectively. Together with the toroidal component, they 
generate a surface of revolution about the polar axis. Taking the 2m(n) power of this poloidal 
flux function, and considering the complement contours, we could generate the mass density 
profiles of p(^) = pofci^) = Po(C — l^l 2 "^™)). The generalized total plasma pressure fol- 
lows with an added p*i = 1/z decay with z. We remark that the nonlinear solution of ©(a;) 
is driven by the generalized total pressure p*. According to this model, although 0(x) is 
occasionally symmetric, this function is basically asymmetric and often with (3 P ^> 1. Conse- 
quently, supernovae are asymmetric in nature with different degree of asymmetry. Observed 
from broadside, this asymmetry should be evident. As the line of sight moves towards the 
poles, the supernova becomes less asymmetrical and more spherical by projection effects. 
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The proto-supernova cavity, filled with circulating plasmas and magnetic fields, generates a 
sequence of rotational equilibria as (3 P increases. When {3 P gets sufficiently high to crack the 
stellar hydrogen envelope, an often asymmetric supernova is erupted resulting in the recoil 
of the neutron star. Since /3 ^> 1, plasma pressure is the prim ary trigger of supernova rath er 
than magnetic field, which is consistent to some simulations (IKomissarov fc Barkovl 120071 ) . 
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5. Collimated Polar GRB Jets 

We can also arrange the Grad-Shafranov equation, Eq. [221 as 



R dr 2 



B at) \smv do J 
= n{n + l) . (28) 



The equation and its solutions are 



;i - x 2 )^M + n{n + l)Q(x) = , (29a) 
ax 

Q(x) = (l-x 2 )^-^- = (l-x 2 ) , (29b) 



where P n (x) is the Legendre polynomials and the last equality is for n — 1. As for the R 
equation, it reads 



r 2 ^rr - n(n + l)R = Tm(3 v z%Az)R 2m - 1 sin 2 6\Q\ 2m - 2 
dr 2 

= TmP^R 2 ™- 1 sin 2 9\0\ 2m - 2 , (30) 

where we have made use of Eq. [TBI to get the last equality. To make this equation separable, 
we consider 2m — 2 = — 1 and then take n = 1 in Eq. I29bl to get 



r 2 — - n(n + l)R = Tmf3 p z 3 R 2 ™- 1 — — = Tm(3 p z s . (31) 
dr 1 B 

With n = 1 and 2m = 1, the solution of R{z) is given by the homogeneous solution in 
bracket plus the particular solution 

R(z) = (az n+l + b^A + cz 3 , (32) 
c = =F^m/3 p . 
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To solve for equilibrium R(z), we again take the lower sign in Eq. [3TJ Taking a = +2, 
b = +2, c = +(3 p /8, and f3 p = 5, Fig. 11 shows the function R(z). Together with 6(x) = 
(1 — x 2 ), the poloidal magnetic field lines and the poloidal density weighed velocity stream 
lines, ty(z,x) = R(z)Q(x) = C, run on the contour line as shown in Fig.12 with C = 
2,3,4,5,6,7. Together with the toroidal component, they generate a surface of revolution 
about the polar axis. The magnetic field lines of this jet in the cavity can be closed via 
the stellar envelope. For large distances, the surface is dominated by the z 3 term with 
ty(z, x) = cz 3 (l — x 2 ) = zz\ = C, where z± is the perpendicular distance from the polar 
axis. Let us approximate the contour as 



C 

Z\\ 



(33) 



This shows that as z increases, z± would decrease, giving a cusp funnel structure, or a 
tornado-like line vortice structure, of a polar jet. In Fig.12, only one jet is shown for more 
clarity, and the opposing jet is likewise. Furthermore, only the initial part of the funnel 
contours are shown. For y{9 = 0) > 15, the analytic funnel solution of Eq. [33] should be 
superimposed on it. We also note that the inner contours 2,3 and 4 are connected to the outer 
magnetosphere, while the outer contours 5,6, and 7 would have to come from external sources, 
like an accretion disk. If these external sources are unable to provide such corresponding field 
strengths of the contours, these outer jet contours should be ignored. For small distances, 
the poloidal flux surface is dominated by the b/z term with ^>(z,x) = b(l — x 2 )/z = C, such 
that z is proportional to sin 2 6. Considering sin 2 9 as the dipole field of the magnetar, the 
polar funnel structure connects up with the magnetosphere of the central compact star, as is 
shown in a close-up look in Fig. 13. We name this as the magnetospheric jet. This connection 
should take place outside the light cylinder of the magnetosphere. 

Furthermore, for a given zn, Eq. [331 shows that z]_ is proportional to the contour value C . 
As a result, the gradient of \I/ increases outward. This means stronger field lines and faster 
stream lines are on the outside. Because of the overall pressure balance, higher pressure 
plasmas are on the inside close to the axis. Plasma confinement is, therefore, accomplished 
by the strong magnetic fields outside surrounding the plasmas. The magnetic confinement is 
also enhanced by the inward magnetic curvature of the cusp funnel geometry which puts a 
magnetic tension force on the plasmas. As for the mass density, which satisfies v- Vp(\I / ) = 0, 



16 



the density gradient Vp(\l/) = (dp/d^f)'V^> is perpendicular to the \l/ contour surface. The 
mass density profiles are given by the complement contours of the poloidal flux function, 



according to p(\I/) = po(C — 



1 2m\ 



Po{C — 1^1) with 2m = 1, and the generalized total 



plasma pressure p* follows accordingly Although the cusp funnel volume is unbounded, the 
jet will be filled up to the point where the energy flux of the star's magnetosphere can supply 
over time. Upon plowing through the stellar envelope along the polar axis, this cavity jet 
could then erupt into a gamma-ray burst jet. 

According to our findings, the outer magnetosphere is torn open by the plasma ram 
pressure, and is flipped polewards by the poloidal plasma flows to form a cusp funnel. Con- 
sequently, the magnetic field lines of this funnel wrap around a high pressure plasma column 
around the axis. Usually, when only toroidal rotational plasma velocity is taken into ac- 
count, the outer magnetosphere opens up into monopole-like radial fields. In our case, the 
poloidal plasma velocity brings the outer magnetosphere to the polar direction . This of- 
fers a MHD description of the magnetic towering mechanism (jLynden-Belll 120031 ) . Standard 
jet formation schemes rel y on the angular momentum of a binary system o r an accretion 



disk-compact star system (IPrendergast 



2005; 



Gourgouliatos fc Vlahakis 



2010|). A single star 



discounting the external stellar envelope can hardly provide the needed angular momentum 
to form jets. The present MHD model can provide asymmetric supernova and polar colli- 
mating gamma-ray burst jet con figurations in the cavity, therefore, favoring the gamma-ray 
burst and supernova association ( 



Soderberg et al. 



2006; 



Galama et al. 1999; 



Mazzali et al. 



Campana et al 



2006 



Pian et al. 



2006 



20061 ). Our axial jet has a specific "beam pattern" . Such 



feature is compatible with models that consider gamma-ra y burst and afterglow properties 



as a r esult of viewing angle on the beam pattern of the jet (ILipunov et al 



2001 



Rossi et al. 



2002 



Salmonson fc Galama 



2002). 
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6. AGN Plasma Torus and Conclusions 

Active galactic nuclei (AGN) is a term used to designate astrophysical objects of Seyferts, 
quasars, radio galaxies, and blazers (BL-Lacs), that emit mighty electromagnetic radiations 
at different bands from radio, optical, to X-ray, gamma-ray. From the spectroscopic char- 
acteristics of these objects, it is deduced that they could be the re sults of a single unified 



AGN structure viewed at different angles. This unified AGN model (jAntonuccilll993l ) would 
consist of a central core with a black hole candidate, an acretion disk, that defines the equa- 
torial plane, and a dust torus further out on this plane. On the rotational (polar) axis of the 
central core, there would be two opposing jets that emit strongly beamed and polarized radio 
waves. Viewed from face-on (down the polar axis) would give blazers, and from edge-on (on 
the equatorial plane) would give type 2 Seyferts, quasars, and radio galaxies with narrow 
emission lines. While viewed at an intermediate angle would give type 1 Seyferts, quasars, 
and radio galaxies with broad emission lines. The material in the acretion disk gravitates 
inward, by transporting its angular momentum outward through dissipations, and feeds the 
central core to generate the jets. The existence of a dust torus in AGN has been inferred 
from high resolution instruments and images. 

To account for this unified AGN structure, we consider another solution of Eq. [32] with 
a negative coefficient a. Writing the negative sign explicitly, we have with n — 1 



R(z) = -az n+1 + [ b— + cz 3 ) . (34) 



z n 



The two terms in bracket generates a minimum in the positive R(z) domain. The first term 
shifts the minimum down to the negative domain bounded by Z\ and Z2 where R(z\) = 
R(z 2 ) = 0. With the negative sign explicit, we take again a = +2, b = +2, c = +(3 p /8, and 
f3 p = 5, the function R(z) is shown in Fig. 14 with three regions, separated by z\ = 1.17 and 
Z2 = 3.09. In the absence of the envelope, the poloidal flux function ^(z,x) = R(z)Q(x) 
with B(x) = (1 — x 2 ) is shown in Fig. 15, which consists of a dipole-like magnetosphere, a 
plasma torus, and a polar jet separated by two spherical separatrix at z\ and z%. This polar 
jet can be formed if there is horizontal accretion to feed it. We name this as the accretion jet 
in contrast to the magnetospheric jet. A close-up look at the AGN magnetosph ere is shown 
in Fig. 16, which gives an intrinsic AGN magnetic moment ( Schild et al. 20061 ) . Naturally, 



upon this basic AGN structure, we can superimpose an AGN magnetospheic jet as discussed 
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in the preceeding section to generate quasars and blazers. 

To conclude, we have used a set of MHD equations with divergence-free axisymmetric 
poloidal and toroidal non-field aligned plasma flows to study the dynamics of the cavity 
between the central core and the external envelope in the final phase of a collapsing star. 
A sequence of steady state rotational MHD equilibria in response to the increasing plasma 
pressure is solved to represent the plasma evolution in the cavity. The spatial configuration 
is described by the rotational Grad-Shafranov equation where the ratio of the generalized 
total plasma pressure to the poloidal magnetic pressure, f3 p = 2fip /(a 2 A ) 2 , is the cavity 
parameter. By assigning two source functions, the rotational Grad-Shafranov equation can 
be solved for asymmetric supernova, polar collimated cusp funnel gamma-ray burst jet, and 
active galactic nucleus plasma torus. It is important to remark that both the asymmetric 
supernova lobes and the cusp gamma-ray burst polar jets are connected directly to the mag- 
netosphere of the central compact star, not to an accretion disk. This structure identifies the 
magnetosphere as the central engine of supernova and gamma-ray burst events, by provid- 
ing plasma, magnetic energy, and rotational energy. Since f3 p gets much larger than unity, 
plasma pressure is likely to be the primary agent in cracking the stellar envelope, instead of 
the magnetic field. 
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x(e=7i/2) 



Fig. 1. — The lobe structures of Q(x) for n — 1, 2, 3 with f3 p = are plotted on the (r — 6) 
plane. 
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Fig. 2. — The 1-lobe structure of Q(x) for n = 2 with j3 p = 9.57 and 0'(— 1, 1) is plotted on 
the (r — 6) plane. 
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Fig. 3. — The 2-lobe structure of Q(x) for n = 3 with j3 p = 6.69 and 0'(— 1, 3) is plotted on 
the (r — 6) plane. 
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Fig. 4. — The 2-lobe structure of Q(x) for n = 3 with f3 p = 6.12 and 0'(— 1, 2) is plotted on 
the (r — 6) plane. 
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Fig. 5. — The 3-lobe structure of Q(x) for n = 4 with (3 P = 6.05 and ©'(—1, 4) is plotted on 
the (r — 6) plane. 
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Fig. 6. — The 2-lobe structure of Q(x) for n = 4 with (3 P = 6.41 and 0'(— 1, 4) is plotted on 
the (r — 6) plane. 
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Fig. 7. — The 3-lobe structure of Q(x) for n = 4 with /3 P = 12.3 and 0'(— 1, 3) is plotted on 
the (r — 6) plane. 
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Fig. 8. — The 3-lobe structure of Q(x) for n = 4 with (3 P = 47.0 and 0'(— 1, 1) is plotted on 
the (r — 6) plane. 
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x(e=7i/2) 

Fig. 9. — The poloidal flux function ty(z,x) contours for n = 3 with f3 p = 6.12 are plotted 
on the (r — 9) plane with contour values of C — 1.5, 1.0, 0.7, 0.5, 0.3 to show the asymmetric 
supernova poloidal magnetic field lines and density weighed plasma velocity stream lines. 



32 




x(e=7i/2) 

Fig. 10. — The poloidal flux function ^(z, x) contours for n = 4 with f3 p = 6.41 are plotted on 
the (r — 6) plane with contour values of C = 4, 3, 2, 1, 0.5 to show the asymmetric supernova 
poloidal magnetic field lines and density weighed plasma velocity stream lines. 




Fig. 11.— The profile of R(z) is plotted with a = +2, b = +2, c = +(3 p /8, and /3 P = 5. 
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Fig. 12. — The poloidal flux function ty(z,x) contours are plotted on the (r — 0) plane with 
contour values of C — 7, 6, 5, 4, 3, 2 to show the colimating GRB cusp funnel along the polar 
axis. 




x(8=7t/2) 

Fig. 13. — The poloidal flux function x) contours are plotted on the (r — 9) plane with 
contour values of C — 7, 6, 5, 4, 3, 2 to show the openning of the magnetosphere to the polar 
cusp. 
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z 



14. — The profile of R(z) with a negative domain is plotted with a = +2, b = +2, 

+/3 p /8, and /3 P = 5. 
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Fig. 15. — The poloidal flux function ty(z,x) contours are plotted on the (r — 9) plane to 
show the AGN magnetosphere, plasma torus with C = —1,-1.5,-2, and polar jet with 
C = 7, 5, 3, 1, separated by two spherical separatrix at z\ = 1.17 and z\ = 3.09. 
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Fig. 16. — The poloidal flux function ty(z,x) contours are plotted on the (r — 9) plane to 
show the AGN magnetosphere with C — 7,5, 3, 1. 



